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RESUMEN 

Favor de proporcionar un resumen en espanol. If you are unable to 
translate your abstract into Spanish, the editors will do it for you. 

ABSTRACT 

We combine archival spectral observations of the LBV star V532 (Romano's 
star) together with the existing photometric data in the B band. Spectroscopic 
data cover 15 years of observations (from 1992 to 2007). We show that the 
object in maximum of brightness behaves as an emission line supergiant while 
in minimum V532 moves along the sequence of late WN stars. In this sence, 
the object behaves similarly to the well-known Luminous Blue Variable (LBV) 
stars AG Car and R127, but is somewhat hotter in the minima. We identify 
about 100 spectral lines in the 3700-^7300A wavelength range. For today, our 
spectroscopy is the most comprehensive for this object. The velocity of the 
wind is derived using Hcl triplet lines (360 ± 30km s _1 ). Physical parameters 
of the nebula around V532 are estimated. 

Key Words: stars: individual: Romano's star (M33) — galaxies: individual: 
M33 — stars: supergiants — stars: Wolf-Rayet 

1. INTRODUCTION 
Luminous Blue Variable s (LBVs) are a class of rare astrophysical objects 



introduced by IContil (| 1984T ) and remain a subject of intense interest. LBVs 
are broadly accepted as very massive and energetic stars emitting close to 
the Eddington limit, evolving from Of toward Wolf-Rayet stars. However, 
contemporary evolutionary models do not answer the question about the exact 
relations between LBV, nitrogen-rich Wolf-Rayet (WN) and hydrogen-rich 
WN (WNH) stars. Numerical simulations show that the objec ts may pass the 



LBV evolutionary stage either before or after the WNH stage ([Smith fc Conti 



2008). Currently, only 35 LBV and LBV candidates are known in our Galaxy 



(jClark et al.l 120051 ) . Studying LBVs in nearby galaxies is very important for 
understanding stellar evolution and the evolution of the interstellar medium 
perturbed and contaminated by massive stars at various evolutionary stages. 

The object V532 (a = 0l' l 35 m 09. s 71, 5 = +30°41'57"1) is located in 
the outer spiral arm of the M33 galaxy. Giuliano Romano was the first to 
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recover its light curve ( Romano! 1 19781) and to find irregular magnitude vari- 
ations between 16™7 and 18™1. Romano classified V532 as a variable of the 
Hubble- Sandage type by the s hape of the light curve and its colour index. 
Humphreys fc Davidson! ( 1994 ) on the basis of its light variations classified 



the star as an LBV candidate. Two maxima of brig htness were detected dur- 
ing the last half of the century (jKurtev et al.l 12001). The first was observed 
around 1970 and the second in the early 1990s (jKurtev et al.l 120011 ). Long- 
term variability has an amplitude of about l m and seems to be superimposed 
on an even stronger downward trend. Short-timescale variability with am- 
plitude ^P 1 ^ was discovered in addition to the longer-timescale variability 
(|Kurtev et al.l 1 2 It ISholukhova et al.l 120021 ) . Such photometric behavior is 



typical for an LBV star. 

First optical spectrum was obtained by T. Szeifert (|Szeifertlll996h at the 
3.5 m Calar Alto telescope in 1992. In this spectrum, "few metal lines are 
visible, although a late B spectral type is most likely (faint Hel)" . This spec- 
trum, obtained with the TWIN spectrograph of the Calar Alto telescope, is 
unique in being obtained in a profound flare state. It is drastically different 
from the hot spectra observed in the minima of brightness (see below) . 

Second spectrum was obtained by O. Sholukhova at the 6m telescope of 
the Special Astrophys ical Observatory (SAO) of Russian Academy of Sci- 
ences (RAS) in 1994 (|Sholukhova et al.lll997j) . Another spectrum was ob- 
tained at SAO with the Multi-Pupil Fiber Spectrograp h (MPFS) in Septem - 
ber 1998. Thi s spect rum was classified as WNlO-WNll lFabrika et al.l |2005h . 



Polcaro et al.l (j2003l ) estimated the bolometric absolute magnitude of the ob- 
ject as Mf,oi ~ — 10'"4, using bolometric correction "of at least -3 mag" and dis- 
tance modulus m — M = 24 ™8. They classify V532 as an LB V because the ob- 
ject fulfills all the criteria of lHum phrevs fe Davidson J 19941 ) . Using five spec- 
tra carried out in 2003 ~ 2006. lviotti et al.l |2006[ l2007t ) find anti-correlation 
between equivalent widths of the Wolf-Rayet blue bump at 4630 ~ 4686A and 
visual luminosity. 

Comparing the spe ctra published by Szeifert ( 1996 ). Fabrika et al. ( 2005! ) 
and lViotti et al.l (|2007l ). we find that the object changes its spectral properties 
significantly. 

Here we combine the spectral observations (both new and already pub- 
lished) with the light curve of the object to trace the spectral variability of 
the star. We describe the data and data reduction process in the next section. 
Results are presented in section [3] and discussed in section 2) In section [5l we 
present the conclusions. 



2. OBSERVATIONS AND DATA REDUCTION 

In this work, we use archival data from the 6m SAO telescope (available 
via ASPID database, |http : //alcor . sao . ru/db/aspid/ 1 and the SUBARU 
telescope, which is operated by the National Astronomical Observatory of 
Japan. The 6m telescope data were obtained with the Multi Pupil Fiber 
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Spectrograph (MPFS) dAfanasiev et alJboOll) and with the SCORPIO multi- 
mode focal reducer in the long-slit mode (jAfanasiev fc Moiseevl 120051 ) . The 
data from SUBARU w ere obtained with the Faint Object Camera (FOCAS) 
(jKashikawa et al.ll20 02) in the Cassegrain focus. Uncertainties are everywhere 
dominated by statistical Poissonian noise (readout noise and round-off errors 
are significantly smaller). 



2.1. MPFS Data 

MPFS obtains simultaneously the spectra from 240 (16 x 15) or 256 (16 x 
16) spatial sampling elements arranged in the form of a rectangular array 
of square lenses. In 2002 and 2004-2005, the data were obtained in the 
240-element and 256-element configurations, respectively. Spatial sampling 
of l"x 1" was used. Light from individual sampling elements is collected by 
microlenses and transmitted by means of optical fibers reorganized in the 
form of a pseudo-slit toward the spectral camera. Grating # 4 (600 Apix -1 ) 
providing spectral resolution of about 6 A was used for all the observations. 
Detectors CCD TK 1024 (1024 x 1024 pixels) and EEV42-40 (2048 x 2048 
pixels) were used in 2002 and 2004-2005, correspondingly. Sky background 
spectrum at the distance of 4' away from the object is taken simultaneously 
with the object by 17 (for 16x16 field size) or 16 (for 16x15 configuration) 
additional fibres. We summarize the relevant information on the MPFS data 
in table [TJ 

Data reduction system was written in IDL environment and makes use of 
procedures written by V. Afanasiev, A. Moiseev and P. Abolmasov. Reduction 
process consists of the standard steps for panoramic data reduction (see for 



example I Sanchez! ((2006)): bias subtraction, flat-ficlding, removal of cosmic-ray 
hits, extraction of the individual spectra from the frames and their wavelength 
calibration using a spectrum of an He-Ne-Ar lamp. At every wavelength, we 
calculate the median sky level using the offset fibres and subtract from the 
spectra of the field. Spectra of spectrophotometry standard stars were used 
for absolute flux calibration. 

Three emission-line spectra obtained with MPFS in 2002-2005 were used. 
We estimate signal-to-noise ratio in continuum as 10-^30 per resolution ele- 
ment (see table [T] for details) . Integral spectra were extracted in an annular 
aperture 2" in radius. Up to the instrumental resolution, the object point- 
like. Note that, unlike long-slit spectrographs, MPFS is free from slit losses 
at absolute flux calibration is much more reliable. 



2.2. SCORPIO Data 

Observational data from SCORPIO are summarized in table [2] All the 
spectra were reduced using ScoRe package for long-slit data reduction, written 
in IDL (available at |http : //narod . ru/disk/5 2380090 007Score_vl . 2 . tar . html[ ). 
CCD frames were de-biased, cosmic particle hits were removed from all types 



4 



MARYEVA & ABOLMASOV 



16 



* 17 



18 



19 



1995 



2000 



2005 



TWIN 



i+ ft 



MPFS 



4 1" 



f 5 + i 




- 4^ f 



SUBARU 



XI in 

4.9x10" 5.0x10" 5.1x10" 5.2x10* 5.3x10 4 5.4x10 4 

JD, 2.4 10 6 + 



Fig. 1. B-band light cur ve of V532. Arrow s indicate when spectroscopic observations 
were made. Data from (|Viotti et alj|2007l ) are shown by diamonds. 



of exposures except bias. Save for the single 10 minute long exposure ob- 
tained in February 2005, several object exposures are present for every day 
that simplifies the removal of cosmic hits. Then we divide object and spectral 
standard exposures over normalized mean flat-field frame. After wavelength 
calibration using He-Ar-Ne lamps and night sky OI (A5577 or A6300 depending 
on grism) emission lin es, the CCD data were flux-calibrated using standard 
stars from lOke I (|l990l ) spectrophotometric standard list (see tabled]). Spectra 
were extracted by fitting the profiles of slices across dispersion with Gaussian 
function. 

In total, eight spectra were obtained with SCORPIO in 2005-2008 (dates 
of observations are shown by upward arrows in figure [1]). Spectral resolution 
is 10 A(for VPHG550G grism) and 5 A(for VPHG1200G and VPHG1200R 
grisms), signal-to- noise ratios in continuum vary from about 10 to about 50 
for larger total exposures. 



2.3. Data from SUBARU 

One exposure 1200s in length was obtained with the SUBARU telescope in 
October 2007. VPHG450 grism was used providing the spectral range of 3750- 
5250 A. Slit width of 0"5 implies spectral resolution of about 1.7 A (for today, 
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TABLE 1 

OBSERVATIONAL LOG FOR MPFS DATA. S/N IS SIGNAL-TO-NOISE 
RATIO PER RESOLUTION ELEMENT. 

Date Exposure Seeing, " S/N Spectral Spectral 

time, s standard star range, A 



05. 10. 2002 900 
13. 11. 2004 4200 
17. 01. 2005 3600 



3.8 16 
1.5 20 
1.5 18 



BD25d4655 4250-6700 
HZ44 4000-7000 
G248 4000-7000 



TABLE 2 

OBSERVATIONAL LOG FOR THE SCORPIO DATA, a IS SEEING, PA 
IS POSITION ANGLE, S/N IS SIGNAL-TO-NOISE RATIO. 



Date 




Exposure 
time, s 


Grism 


Spectral 
range, A 


SX, 
A 


S/N 


a, " 


Spectral 
standard 
star 


PA,° 


6. 02. 


2005 


600 








8 


1.7 


G248 


-136 


30. 08 


. 2005 


1200 


VPHG550G 


3500-7200 


10 


30 


1.9 


G191-B2B 


210 


8. 11. 


2005 


3300 








45 


1.9 


BD25d4655 


145 


3. 08. 


2006 


1500 








20 


2 


GD248 


200 


10. 08 


. 2007 


1800 


VPHG1200G 


4000-5700 


5 


24 


2 


BD33d2642 


252 


5. 10. 


2007 


2700 








40 


1.1 


BD25d4655 


-141 


8. 01. 


2008 


1800 


VPHG1200R 


5700-7500 


5 


20 


2.1 


BD25d4655 


48 


10. 01 


. 2008 


1800 


VPHG1200G 


4000-5700 


5 


20 


1.4 


BD25d4655 


18 
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it is the best resolution achieved for this object) in the extracted spectrum. 
Data were reduced using IDL-based software. The CCD frames were bias- 
subtracted and fiat-fielded. We used the Th-Ar arc spectrum for wavelength 
calibration. Star BD40d4032 from SUBARU spectrophotometric standard list 
was used to calibrate the stellar spectra. We extract the spectrum in the same 
way we did it for SCORPIO (see above section [2~2j) 

If we degrade the spectral resolution to 5A, the spectrum becomes prac- 
tically identical to the spectrum obtained with SCORPIO on the same date. 
Though the spectral shapes are similar, the spectra differ in normalization 
(about a factor of three) connected to slit losses. 



2.4. Photometric Data 



The light curve (see figured]) was provided by Vitalij Goranskij and consists 
of the CCD data obtained by Goranskij and Zharova with the lm SAO tele- 
scope and two instruments of the Crimean laboratory of the Sternberg Astro- 
nomical Institute (SAI) and photographic plates from the SAI collection. CCD 
data were reduced with Maxim DL softwar e (http : //www . cyano gen . com/maxim_main . php) 
The joint light curve will be published by Zh arova et al. ( 201dh in a separate 
paper, all the details of the data reduction process will be given there. Photo- 
graphic B-band magnitudes were identified with the B-band magnitudes ob- 
tained by CCD observations. CCD data uncertainties are of the order 0™05, 
plate data have larger errors depending on the source brightness, usually of 
the order 0.1 - 0.2 m . 



The joint curve contains the data of IViotti et al.1 (|2007h but complements 
it with a much larger photometric material and primarily allows to trace the 
recent evolution of the object. 



3. RESULTS 

3.1. Spectral Evolution 

The first optical spectrum of V532 was obtained during the rise of bright- 
ness, when the object was 16™4 in the V band. This spectrum covers two 
spectral ranges AA4400 — 5000 and AA5800 — 6750, as it was obtained with the 
two-channel TWIN spectrograph (a comprehensive description is present at 



http: //www. caha.es/CAHA/Instruments/TWIN/HTML/twin.html). In this 



spectrum, two strong lines are present, Ha (shown in figure [2]) and Hf3, 
having complex profiles consisting of a narrow line with P Cyg profile and 
broad wings. FWHM (Full Width at Half Maximum intensity) of the wings 
is ~ 15 A for H/3 and ~ 20 A for Ha. Such wings were fir st found for LBVs 



in the spectrum of P Cygni by iBernat fc Lambertl (|1978l ) . These wings are 



explained by scattering of line photons by free electrons in the stellar wind. 
Similar line profiles were o bserved for the LBV star s R127 and AG Car during 



the initial rise to maxima ( Stahl et al.l ll983, 2001) 
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TABLE 3 



LIST OF EMISSION LINES DETECTED IN THE FOCAS AND 
SCORPIO (JANUARY 2008) SPECTRA OF V532. FOR HIGHER 
SIGNAL-TO-NOISE RATIOS, EQUIVALENT WIDTHS (EW) ARE 
GIVEN. FOR LINES WITH P CYG PROFILES, WE GIVE BOTH 
EMISSION AND ABSORPTION COMPONENT EWS. 



A, A 


Ion 


EW 
emission, 
A 


EW 
absorption, 
A 


A, A 




Ion 


EW 
emission 
A 


3770.60 


Hll+Hell 






4613 


.90 


Nil 




3797.90 


HlO+Hell 






4621 


.40 


Nil 




3819.76 


Hel 






4630 


,54 


Nil 




3835.39 


H9+HeII 






4634 


.00 


NIII 


6.2 ±1.0 


3871.82 


Hel 






4640 


.64 


NIII 


7.4 ±0.5 


3889.05 


H8+Hcl 


11 ± 1 


3.7 ± 1.3 


4643 


.09 


Nil 




3964.73 


Hel 






4650 


.16 


C III 


5.0 ±0.5 


3970.08 


He+HcII 






4658 


.10 


[Felll] 


2.5 ± 1.0 


3994.99 


Nil 






4685 


.81 


Hell 


15.2 ±0.2 


4009.00 


Hel 






4701 


.50 


[FcIII] 


2.9 ± 1.0 


4025.60 


Hcl+Hell 


2.5±0.5 


1.2±0.2 


4713 


.26 


Hel 


1.6 ±0.2 


4088.90 


SilV 






4861 


.33 


H/3 


22.5 ± 1.5 


4097.31 


N III 






4921 


.93 


Hel 


4.3±0.2 


4101.74 


H(5+HeII 


7.6±1.0 




4958 


.91 


[O III] 


1.3±0.2 


4103.40 


NIB 






5006 


.84 


[O III] 


3.9±0.6 


4116.10 


SilV 






5015 


.67 


Hel 


2.6±0.3 


4120.99 


Hel 






5411 


.50 


Hell 


6.2±0.2 


4143.76 


Hel 






5666 


.60 


Nil 




4199.80 


Hell 






5676 


.02 


Nil 




4236.93 


Nil 






5679 


.56 


Nil 




4241.79 


Nil 






5686 


.21 


Nil 




4241.79 


Nil 






5710 


.76 


Nil 




4340.47 


H 7 +HeII 


8.0±0.3 




5875 


.79 


Hel 


25.5±0.7 


4387.93 


Hel 


1.4±1.0 




6548 


.00 


[N II] 




4471.69 


Hel 


5.0±0.3 


1.7±0.2 


6562 


.82 


Ha+Hell 


107±3 


4481.13 


Mgll 






6583 


.00 


[N II] 




4510.90 


NIII 






6678 


.15 


Hel 


18 ±3 


4514.90 


NIII 






6683 


.20 


Hell 




4518.20 


NIII 






7065 


.44 


Hel 


18 ±3 


4523.60 


NIII 






7135 


.73 


Arlll 


3.2 ±0.3 


4530.80 


NIII 






7281 


.35 


Hel 




4534.60 


NIII 














4541.60 


Hell 














4547.30 


NIII 














4601.50 


Nil 














4607.20 


Nil 
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Fig. 2. Ha line profile in the Calar Alto spectrum obtained in 1992. The spectrum 
is normalized by the local continuum level. 
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TABLE 4 

SPECTRAL CLASSES IDENTIFIED VIA THE SCHEME OF SMITH, 
CROWTHER AND PRINJA (1994) 



B, mag Spectral Date B, mag Spectral Date 





subtype 






subtype 




17.5 


WN10.5 


2002/10/05 


17.6 


WN9 


2005/11/08 


16.9 


WN11 


2004/11/13 


18.3 


WN8 


2006/08/03 


17.1 


WN11 


2005/01/17 


18.4 


WN8 


2007/08/10 


17.15 


WN11 


2005/02/06 


18.5 


WN8 


2007/10/05-08 


17.3 


WN10 


2005/08/30 




WN8 


2008/01/08-10 



Emission lines of HcIA5876 and SiI IA5957.612, 59 78.970, 6347.091, 6371.359 
are also present in the red spectrum. Szeifert ( 19961 ) mentions these Sill lines 
as "weak metal emissions" . 

Figure [3] shows all the spectra of V532 in the blue range (4000-5500A) 
analysed in this article, obtained at different spectral resolutions. All the 
spectra are normalised to continuum level in a uniform way. In order to get 
a reasonable normalisation for our spectra, we chose several wavelength in- 
tervals practically free from Wolf-Rayet emissions (4250-^4270, 5100-^5200, 
5520-^5620, 5750-^5800 and 6950^6970A) and reconstructed the continuum 
using a second-order polynomial. Below, for spectral classification we use 
characteristic equivalent width ratios, that makes our results practically inde- 
pendent of spectral resolution. 

The spectral appearance of the spectra of V532 obtained in 2002-2008 re- 
sembles that of late WN stars. It allows to apply a classification scheme used 
for WN stars, bearing in mind that abundances of individual elements may dif- 
fer, but physical conditions are similar. All the spe ctra obtained between 2002 
and 2 008 were classified using the classification of ISmith. Crowther fc Prinia 
(|1994l ) for WN6-11 stars based primarily on relative strengths of NV AA4604 — 
20, NIV A4058, NIII AA4634 - 41 and Nil A3995 emission lines. The method 
has low dependence on elemental abundances, though only helium and nitro- 
gen lines (preferably, ratios of the lines of one element) are used. The results 
of spectral classification are given in table |U 

Light curve exhibits a local maximum in 2004 and early 2005. In the 
spectra obtained in this period, one may observe strong emission lines of 
hydrogen and neutral helium. The spectral appearance of V532 shows strong 
similarities with a spectrum of a WN11 star. Hel A4713 line is stronger 
than the Nil blend, and HcIIA4686 is absent. In the spectrum obtained in 
2002, the Nil A3995 line is outside the spectral range, therefore we carry out 
classification comparing NIII AA4634 — 41 and Nil AA4601 — 43 blends having 
similar intensity to the Nil A3995 emission. NIII AA4634 — 41 lines appear, but 
are weaker than Nil AA4601 — 43. HeIIA4686 line is approximately as bright 
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4000 



4500 



5000 



5500 



Fig. 3. Optical spectrum evolution in the blue spectral range (4000^52001). Spectra 
are normalized by the local continuum level and vertically shifted for clarity. 



as the Nil AA4634 — 41 emission, hence we classify the object as WN10.5. 

Starting from the middle of 2005, Romano's star weakens in the opti- 
cal range, its visible magnitude falls from 17 to 18™8. During this period 
its spectrum evolves from WN10 (August 2005) through WN9 (November 
2005) to WN8 (August 2007) by three spectral sub-classes. In October 2007, 
Hell A4686 line is stronger than NIII AA4634 - 41. However, NIV lines do not 
appear in the spectrum. We also classify the spectrum as WN8. 

V532 evolution on the equivalent width diagram of Hcl A5876 versus 
Hell A4686 is shown in the left panel of figure |U Locations of well-proven 
Galactic and Large Magellanic Cloud (LMC) Wolf-Rayet stars R84 (WN9), 
Sk-60 40 (WN10), LBV HDE 269582 (WN10) and some other objects are 
shown for comparison with our object. The figure shows that in the period 
between 2002 and 2005, when the star brightened by about one magnitude, its 
spectrum changed from WN9.5 to WN11. Since 200 5, the spectrum changes 
smoot hly according to the sequence established by ISmith. Shara fc Moffat 
dl996h . 

We used a quantitative chemistry-independent criterion based on the FWHM 
of the Hell A4686 line for alternative spectral classification. In Fig01 we show 
the location of V532 on the diagram of the equivalent width of Hell A4686 
versus FWHM of this line. We have measured equivalent width and FWHM of 
the Hell line in the FOCAS spectrum. For this, we approximate Wolf-Rayet 
blue bump by 7 Gaussians. The position of V532 is fully consistent with its 
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Fig. 4. Left panel: V532 location on equivalent width diagram of HeIA5876 versus 
HeIIA4686 for different dates of observation. The points are connected with dotted 
lines in chronological order. Right panel: V532 location in October 2007 (FOCAS 
data) on equivalent width versus FWHM diagram for the HeIIA4686 line. In both 
graphs, known Galactic (open symbols) and LMC (filled) WN stars are shown for 
comparison: WN8 by circles, WN9 by squares, and WN10 by triangles. Data on 
these objects were taken from lSmith et al. (1995). 



Galactic and LMC WN9 analogues. Spectral class defined from the diagram 
is consistent with that determined from the relative strengths of Nil, NIII and 
Hel with accuracy about one subclass. 

Table [3] shows the lines detected in the spectra obtained with FOCAS 
(October 2007) and with SCORPIO (January 2008, with grisms 1200G and 
1200R). Three spectra are used in order to cover the maximal wavelength 
range at maximal possible spectral resolution. The spectrum of V532 does 
not change noticeably from October 2007 to January 2008. Equivalent widths 
of the principal lines (absorptions, emissions and P Cyg) are given with errors 
of approximation. Uncertainties due to the choice of the continuum level are 
smaller than the errors of approximation. Lines with P Cygni profiles were 
resolved by fitting with the models described in the next subsection. 

3.2. Line Profiles and Terminal Wind Velocity 

In figure [5] we show the FOCAS spectrum in the optical blue range 3800 — 
5100A. The Wolf-Rayet blue bump (consisting primarily of NIIIAA4634, 4640, 
CIIIA4650, HeIIA4686 and HeIA4713) is clearly seen in this spectrum. We 
suppose that the lines near 4658 and 4701 A are nebular lines of [Felll]. They 
can not belong to CIV since CIV AA5812, 5801 lines are not present in the 
SCORPIO spectrum obtained simultaneously in October 2007. 

Analyzing the FOCAS spectrum of V532, we found that triplet and sin- 
glet lines of Hel have different profiles (see figure |6]). Triplet lines of Hel 
(A3889, 4025, 4471) show strong P Cyg profiles, while singlet lines (A3965, 4922, 5016) 
have flat-topped profiles. Widths of these lines correspond to velocity span of 
about 100 km s _1 . 
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Fig. 5. Optical spectrum of V532 obtained with FOCAS in the range 3800-5100AA. 
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We used P Cyg profiles of triplet Hel lines to estimate the terminal wind 
velocity Voo . Line profiles were fitted with a sum of two Gaussians (one repre- 
senting the absorption, the other the emission component). We suppose that 
the widths of emission and absorption components are equal to the instru- 
mental profile width of 1.8A which is almost independent of wavelength. The 
instrumental profile width is determined using wavelength calibration spec- 
trum lines of similar wavelengths. Terminal wind velocity is estimated by the 
velocity shift between the Gaussian centers. Even for the emission component 
of Hel lines, we do not detect any Doppler broadening. This may mean that 
the profiles are not true P Cyg but wind blucshiftcd absorptions plus nebular 
emissions. This possibility does not however alter the wind velocity estimates. 

Wind velocities for all the Hel lines with P Cyg profiles are equal within 
the statistical errors. The mean wind velocity for the three triplet lines is 
360 ± 30 km s" 1 . This value is consistent with the t erminal wind velocities 
for late WN (jCrowther et al.lll995t ISmith et al1ll995l ) that we give in table [5] 
for comparison. 



It is interesting to compare these profiles with the HeIIA5412 line in the 
SCORPIO spectrum obtained in October 2007. The spectrum has very high 
signal-to-noise ratio that allows reasonable terminal wind velocity estimates 
in spite of the significantly worse spectral resolution. The line has a classi- 
cal low-optical-depth P Cyg profile: broad emission and blueshifted narrow 
absorption. We approximate this line by a classical P Cyg model profile (flat- 
topped emission with a narrow blucshiftcd absorption component) convolved 
with instrumental point spread function that we consider Gaussian (the in- 
strumental profile width is ~ 5.5A). Similar line profiles are produced by 
optically thin envelopes expanding at a constant velocity. Its fitting yields 
the wind speed of 200 ± 15 km s _1 . The expansion velocity indicates that the 
line is formed in hotter inner parts of the wind that move with lower outward 
velocities. Among Pickering emissions, we also detect A4541.6 and A4199.8 
having similar profiles in the FOCAS spectrum, but we lack signal-to-noise to 
obtain reliable expansion velocity estimated for these lines. 



3.3. Nebular Lines 

The lines at 4959 and 5007 A are clearly identified as [OIII] emissions, 
that is consistent with their flux ratio equal to 3 and FWHMs equal withi n 



the measurement errors. In the cooler spectra of 2003 (jPolcaro et al.ll2003h . 
nebular lines of [OIII] are overridden by Nil 4994-5005 emissions. However, 
in 2007, Nil lines are absent, but the [OIII] A4959, 5007 doublet is clearly seen. 

Nebular lines of [OIII] AA4959, 5007, [Nil] AA6548 - 83 are present in all 
the spectra analysed by us while [SII] A6717 — 31 doublet is not. Therefore, we 
can estimate the electron density of the surrounding nebula. Electron density 
is between 10 3 6 cm -3 and 10 4 9 cm~ 3 . The former value corresponds to the 
critical density of the [SII] doublet, the second to that of [Nil] AA6548 — 83 



lines (see for example lOsterbrock fc FerlandJ (|2006f )) 
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Fig. 6. Top panel: profiles of triplet HeIA3889, 4025, 4471 lines (from left to right). 
Two-component model fits are shown (see text for details). Bottom panel: profiles 
of singlet HeIA3965, 4922, 5016 lines. FOCAS data. 

All the nebular lines present in the FOCAS spectrum as well as the 
[Arlll] A7135.73 emission in SCORPIO data show flat-topped or two-peaked 
profiles. The widths of emission line cores are of the order 100 km s _1 . Once 
we know the density of the emitting gas, we may estimate the size of the 
Stromgren region and the mass of the ne bula surroun ding V532. The former 
may be found as follows (see for example iLena (|l974l )): 



R s 



\4iran2 J 



1/3 



2.6 • 



where a is Case B recombination coefficient for hydrogen, a 

/ 4 \ 0.85 

10 -13 I ) cm 3 s _1 , n e is electron density, and S is the number of hydrogen- 
ionizing quanta production rate. For 09.51a stars, prob ably similar to the ob- 
ject i n mass and luminosity, S is estimated as 10 49 ' 17 s _1 (jOsterbrock fc Ferland 
2006T ). One may estimate the size of a homogeneous nebula around V532 as: 



R s -0.1 



( s 



\ 10 49 s 



1/3 



( n * ) 

V10 4 cm- 3 / 



-2/3 / t y 
V 10 4 K ) 



pc 



(1) 



Note that possible inhomogeneity has little effect on the linear size: for a 
gas filling only some part of the nebular volume /, R s should scale with filling 
factor as cx / _1 ^ 3 . However, sizes of about lpc are plausible. Nebular mass 
may be trivially inferred as: 



47T - 



0.17 ■ 



( s 



-) ( ne r 1 (—) 

J Vl0 4 cm- 3 / V lo4K / 



0.84 



Alp 



(2) 



V lO^s- 1 J V 10 4 cm" 3 J V 1Q4K / 
Estimated size and the ionised mass of the nebula arc consistent with these 



for ejecta of LBV stars (see lSmith. Crowther fc Prinial (|1994l ) and references 
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therein) by order of magnitude. It is possible that the total mass lost by the 
object is higher, about several Solar masses, but we observe only the ionized 
part of the envelope. The emitting gas was probably ejected during one or 
several outburst events at wind velocities about 100km s _1 . 



4. DISCUSSION 



In tabled terminal wind velocities for the WN9-WN11 stars in the Galaxy, 
LMC and M33 and LBV star AG Car (during minimum) are given. These 
velocities are also estimated through the opt ical Hel lines. Because wind 
accel eration is tightly connected to metallicity (jMassev et al I Puis et al.l 

20081 ). we also give oxygen abundances (12 + lg O/H) for selected HII regions 



adjacent to MCA1-B, B517 and BE381. They may be used as a measure for 
initial object mctallicities. 

V532 is located at a distance of about 17' from the centre of M33. Taking 
the central oxygen abundance value of 12 + log O/H = 8.36 ± 0-04 a nd radial 
gradient of -0.027 ± 0.012 dex kpc" 1 (jRosolowskv k SimonjEooj ). we may 



estimate the primordial oxygen abundance for V532 as 8.26 ± 0.08. Ambient 
metallicity is thus similar to that for the 30 Dor star-forming region. Our 
estimate for the terminal wind velocity of V532 is fairly consistent with this 
for other late WN stars in similar environment. 

Wolf-Rayet stars are believed to be more evolved objects than Ofpe and 
Be supergiants. Information on elementary abundances may cast some light 
upon this issue. However, to recover reliable abundance estimates, one needs 
sophisticated modelling of both the moving atmosphere of the star and its 
structure. Here, we restrict ourselves only to some semi-qualitative estimates 
for the H/He abundance ratio. We compare the equivalent width ratios of H/3 
to he lium lines in the spectr a obtained in 2007 to those for several WN9 (taken 
from lCrowther et al.l ([1995)) to estimate the relative abundances of hydrogen 
and helium. The equivalent width ratios are given in table |6] Table shows 
that they are similar to those for BE381 (Brey 64 ) that has H /He = 2 and 
is classified as WN9h bv lCrowther k Smith! |l996l) . By analogy with BE381, 
we suppose that for V532, H/He ~ 2. It places Romano's star in the region 
occupied by hydrogen-rich late WN stars that have atmospheres moderately 
contaminated by helium. 



We have already mentioned the strikingly different behaviour of triplet 
lines of Hel having P Cyg-like profiles and singlet emissions. Probably singlet 
lines have weak absorptions that we can not detect due to insufficient signal- 
to-noise or spectral resolution. Besides this, HcIA5015 is contaminated by the 
oxygen [OIIIJA5007 emission. 

Similarity of the line profiles of singlet Hel, [OIII] AA4959, 5007 and [Arlll] A7136 
suggests that both singlet Hel and forbidden emissions are produced mainly 
in the low-density ejecta, expanding at velocities of about 100km s _1 . Prob- 
ably, singlet and triplet lines of neutral helium are formed in different places. 
Because the lowest possible level for triplet lines is metastable, concentrations 



16 



MARYEVA & ABOLMASOV 



TABLE 5 

WIND VELOCITIES ESTIMATED THROUGH OPTICAL HEI LINES 
AND AMBIENT OXYGEN ABUNDANCES. 

Star Galaxy WN Voo 12 + logO/H HII 

subtype km s _1 abundance region 



V532 


M33 


9 


360 


8.26 ±0.08 


(during minimum 










epoch 2007) 










MCA1-B a 


M33 


9 


420 


8.315 ±0.061 


B517 b 


M33 


11 


275 


8.334 ±0.083 


Sk-66°40 a 


LMC 


10 


300 




R84 a 


LMC 


9 


400 




BE381 C 


LMC 


9 


280 


8.37 


WR105 a 


MW 


9 


700 




AG Car 










(during minimum 


MW 


11 


300 




epoch 1985-1990) d 










AG Car 










(during minimum 


MW 


11 


200 




epoch 2002) d 










Smith et all ( 19951) 



30 Dor* 



Crowther et alj (1993) 



jCrowther etafl jl995h 
feroh et al.1 (2009h 
Rosolowskv fc SimonJ {2008) 



Rosa fc MathisJ (|l987h 



TABLE 6 

EQUIVALENT WIDTHS OF HELIUM LINES (IN H/3 UNITS) AND 
HELIUM-TO-HYDROGEN ABUNDANCE RATIOS FOR TWO LATE 
HYDROGEN-RICH WNS AND V532 (2007). 



HeI±H8 H7 Hell Hel Hel 











3889 




4686 


4713 


5876 


R84 


WN9 


H/Hc= 


2.5 


3 


2.3 


3.5 


12.6 




BE381 


WN9h 


H/Hc= 


2 


2.1 


2.4 


1.9 


16.6 


0.8 


V532 








2 


2.75 


1.4 


13.4 


0.8 
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of atoms at the lower levels of triplet transitions (at therefore the intensities 
of absorption line components) should be higher. 

We classify the observed evolution of the object as S Dor variability cycle. 
The main difference with AG Car and other LBVs is that in the optical mini- 
mum V532 becomes much hotter and may be classified as a WN9/WN8 star, 
while AG Car stops at about WN11. Non- monotonous spectral changes and 
variability timescales indicate that these spectral changes hardly correspond 
to any real stellar evolution, rather being connected to the ordinary LBV vari- 
ability. It also means that LBV stars may have spectral classes as early as 
WN8. To our notion, there is only one example of an object that acquired 
even earlier spectral classes of WN6/7 during LBV variability cycle, HD5980 
fearba et al.lll997i: uXoenigsberger fc Morenolboosl ) in SMC. The hot spectra 
of this object are probably connected to lower ambient metallicity, too. 

V532 is one more example of an LBV star temporarily becoming a WN. We 
predict that, vice versa, some of the Galactic or extra-galactic Wolf-Rayet stars 
may prove to be dormant LBVs. Probably, the hottest possible temperature 
for an LBV decreases with metallicity, but more data on the brightest stars 
in different environments are needed. 



5. CONCLUSIONS 



Our results show that the object changes from a B emission line super- 
giant in the optical maximum, through Ofpe/WN (WN10,WN11) to WN9 
and further t oward s a WN8 star in deep minimum. We confirm the result of 



Viotti et al.l (120071 ) that there is an anti-correlation between the visual lumi- 



nosity and the temperature of the star, at larger amount of spectral data. 

V532 spans a large range of spectral classes, becoming one of the first 
(probably, the second, after HD5980) LBV stars noticed to make an excursion 
as deep as WN8 into the Wolf-Rayet domain. It is interesting to trace further 
the evolution of the object. Further monitoring of V532 as well as other late 
WN stars is necessary to establish the evolutionary link between WRs and 
LBVs. Some WNs may prove to be dormant LBVs. 



We are grateful to Vitalij Goranskij, Elena Barsukova and Alia Zharova for 
providing us with photometric data and Olga Sholukhova and Thomas Szeifert 
for the spectrum obtained with TWIN Calar Alto spectrograph in 1992. We 
would also like to thank the anonymous referee for valuable comments and 
drawing our attention to HD5980. 
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